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ABSTRACT 

Context. 

Aims. We present spatially resolved high-spectral resolution A"-band observations of the red supergiant Betelgeuse (a Ori) using 
AMBER at the Very Large Telescope Interferometer (VLTI). Our aim is to probe inhomogeneous structures in the dynamical atmo- 
sphere of Betelgeuse. 

Methods. Betelgeuse was observed in the wavelength range between 2.28 and 2.31 yum with VLTI/AMBER using baselines of 16, 32, 
and 48 m. The spectral resolutions of 4800-12000 allow us to study inhomogeneities seen in the individual CO first overtone lines. 
Results. Spectrally dispersed interferograms have been successfully obtained in the second, third, and fifth lobes, which represents 
the highest spatial resolution (9 mas) achieved for Betelgeuse. This corresponds to five resolution elements over its stellar disk. The 
AMBER visibilities and closure phases in the A"-band continuum can be reasonably fitted by a uniform disk with a diameter of 
43.19 ± 0.03 mas or a limb-darkening disk with 43.56 ± 0.06 mas and a limb-darkening parameter of (1.2 ± 0.07) x 10~'. These 
AMBER data and the previous AT-band interferometric data taken at various epochs suggest that Betelgeuse seen in the K-bmd con- 
tinuum shows much smaller deviations from the above uniform disk or limb-darkened disk than predicted by recent 3-D convection 
simulations for red supergiants. On the other hand, our AMBER data in the CO lines reveal salient inhomogeneous structures. The 
visibilities and phases (closure phases as well as differential phases representing asymmetry in lines with respect to the continuum) 
measured within the CO lines show that the blue and red wings originate in spatially distinct regions over the stellar disk, indicating an 
inhomogeneous velocity field which makes the star appear different in the blue and red wings. Our AMBER data in the CO lines can 
be roughly explained by a simple model, in which a patch of CO gas is moving outward or inward at velocities of 10-15 km s _1 , while 
the CO gas in the remaining region in the atmosphere is moving in the opposite direction at the same velocities. Also, the AMBER 
data are consistent with the presence of warm molecular layers (so-called MOLsphere) extending to ~ 1.4—1.5 R+ with a CO column 
density of ~ 1 x 10 20 cnT 2 . 

Conclusions. Our AMBER observations of Betelgeuse are the first spatially resolved study of the so-called macroturbulence in a 
stellar atmosphere (photosphere and possibly MOLsphere as well) other than the Sun. The spatially resolved CO gas motion is likely 
to be related to convective motion in the upper atmosphere or intermittent mass ejections in clumps or arcs. 

Key words, infrared: stars - techniques: interferometric - stars: supergiants - stars: late-type - stars: atmospheres - stars: individual: 
Betelgeuse 



1. Introduction 

Red supergiants (RSGs) experience slow, intensive mass loss up 
to 1 0~ 4 Mq yr~ 1 . Despite its importance not only in stellar evolu- 
tion but also in the chemical enrichment of the interstellar mat- 
ter, the mass loss mechanism in RSGs is not well understood. 
While radiation pressure on dust grains is often considered to be 
the driving mechanism of mass loss in cool evolved stars, it is 
not clear where and how dust forms in RSGs and how mass out- 
flows are initiated. Alternative scenarios include Alfven-wave- 

Send offprint requests to: K. Ohnaka 

* Based on AMBER and VINCI observations made with the Very 
Large Telescope Interferometer of the European Southern Observatory. 
Program IDs: 080.D-0236 (AMBER Guaranteed Time Observation), 
60.A-9054A, and 60.A-9222A. 



driven winds (Airapetian et al. 120001 Schroder & Cuntz 2005 
120071 1. combination of Alfven waves and the wave damping due 
to dust (Vidotto & Janteco-Pereira 2006), and convective turbu- 
lence combined with radiation pressure on molecules (Josselin 
& PlezH5D7}. 

The atmosphere of RSGs exhibits complicated structures. In 
the lower photosphere, vigorous convective motion is expected 
with the convective cell size possibly comparable to the stellar 
radius (Schwarzschild [T975l Freytag et al. 2002 ). The photomet- 
ric variabilities, as well as the variations in the radial velocities 
of the metal lines in the visible, can be interpreted as caused 
by such giant convective cells (Kiss et al. 120061 Gray [2008). 
In the outer region, extended chromospheres exist. For exam- 
ple, the UV observations of the M supergiant Betelgeuse (a Ori) 
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with the Hubble Space Telescope revealed that the hot (~6000- 
8000 K) chromospheric plasma is more than twice as extended 
as the photosphere measured in the near-IR with a bright fea- 
ture (Gilliland & Dupree [19%! Uitenbroek et al. [1998b . The 
image in the Ha line is even more extended, 4-5 times as large 
as the photosphere (Hebden et al. 1987), consistent with the ex- 
tended chromosphere. However, radio continuum observations 
of Betelgeuse with the Very Large Array (VLA) suggest that 
much cooler (~ 1000-3000 K) gas extends to several stellar radii, 
showing an irregular structure (Lim et al. [1998 ). Non-spherical 
shape of the outer atmosphere of Betelgeuse was also detected 
by mid-IR interferometric observations by Tatebe et al. (|2007 ). 
Furthermore, the narrow-slit spectroscopy of Betelgeuse in the 
10 jum region by Verhoelst et al. (2006) revealed that silicate 
dust forms only at large distances from the star (^20 and 
that AI2O3 may form as close as ~2 This means that the hot 
chromospheric plasma, cooler gas, and AI2O3 dust may coexist 
within several stellar radii from the star, but the cooler compo- 
nent is much more abundant compared to the chromospheric gas, 
because it dominates the radio emission. 

The presence of cool gas in the outer atmosphere of RSGs 
is consistent with dense molecular layers existing close to the 
star, the so-called "MOLsphere", which was proposed by Tsuji 
d2000al |2000b ) to explain the IR spectra of the early M super- 
giants a Ori and yU Cep. While these stars were deemed to be too 
hot for H2O to form, he showed that dense H2O gas with column 
densities of the order of 10 20 cirT 2 and temperatures of 1500- 
2000 K at ~ 1.3-2.0 R+ can explain the spectral features at 2.7 
and 6 /mi which cannot be reproduced by non-gray hydrostatic 
photospheric models. Near- and mid-IR interferometric stud- 
ies also lend support to the presence of the MOLsphere toward 
RSGs (e.g., Pen-in et al. l2004ll2005ll2007t Ohnaka|2004il Tsuji 
2006), although the current, crude MOLsphere models cannot 
reproduce the H2O absorption lines observed at 12 /jm (Ryde et 
al. l2006aH2006bl l. On the other hand, Verhoelst et al. ( 120091 ) ar- 
gue against the presence of such dense molecular gas in the outer 
atmosphere of RSGs. They instead propose that the 2.7 and 6 //m 
features in RSGs can be explained by a continuous (i.e. feature- 
less) dust opacity source such as amorphous carbon and metallic 
iron, although it is not clear whether or not such grain species 
indeed form in oxygen-rich environments. 

A better understanding of the inhomogeneous structure of 
the outer atmosphere of RSGs is a key to unraveling the mass- 
loss mechanism in these stars. Inhomogeneities over the stel- 
lar surface were detected by high spatial resolution imaging of 
a few nearby RSGs. The high-resolution images of Betelgeuse 
at 0.7-1.25 fim with spatial resolutions of down to 30 mas 
(the stellar angular size is ~50 mas at these wavelengths) show 
the wavelength-dependent appearance of asymmetric structures 
(Burns et al. [19971 Tuthill et al. [19971 Young et al. [2000b . 
However, their origin is by no means clear. They may be related 
to large convective cells predicted to be present in cool lumi- 
nous stars (Schwarzschild |T975l Freytag et al. 2002) or alterna- 
tively to thermal instability taking place in the outer atmosphere. 
For example, the magnetohydrodynamical simulation of Suzuki 
( 120071 ) for red giant branch (RGB) stars, which are much less lu- 
minous than RSGs, shows that thermal instability leads to "struc- 
tured" stellar winds with many bubbles of hot gas (~ 10 5 K) em- 
bedded in cool winds (~l-5xl0 3 K). For cooler RSGs, the for- 
mation of molecules may also promote such thermal instability, 
and in particular, CO is an important coolant in the atmosphere 
of late-type stars (Cuntz & Muchmore 1994). 

To glean clues to the origin of the inhomogeneities and the 
mass-loss mechanism in RSGs, high spatial resolution observa- 



tions in IR molecular lines are very effective. The high spectral 
resolution (/l/A/l= 12000) of the near-IR interferometric instru- 
ment AMBER (Astronomical Multi-BEam combineR) at VLTI 
allows us to resolve the CO first overtone lines and to spatially 
resolve inhomogeneous structures within each CO line. In this 
paper, we present high-spectral and high-spatial resolution K- 
band AMBER observations of the prototypical RSG Betelgeuse 
(Ml-2Ia-Ibe). 

2. Observations 

2.1. AMBER 

AMBER (Petrov et al. l2007l > operates in the J, H, and K bands 
with spectral resolutions of 35, 1500, and 12000, combining 
three 8.2 m Unit Telescopes (UTs) or 1 .8 m Auxiliary Telescopes 
(ATs). AMBER is a spectro-interferometric instrument which 
records spectrally dispersed fringes on the detector. With the 
maximum baseline length of 130 m currently available at VLTI, 
spatially resolved spectroscopy with an angular resolution of 
down to 2 mas is possible with AMBER, which enables us to 
study the wavelength dependence of the size and shape of the 
object. AMBER observations with three telescopes allow us to 
measure three visibilities and three differential phases (DPs), as 
well as one closure phase (CP). Visibility is the amplitude of 
the Fourier transform (complex function) of the object's inten- 
sity distribution in the plane of the sky and contains information 
about the size and shape of the object. On the other hand, the 
phase of the Fourier transform (also called "Fourier phase" to 
avoid confusion) contains information about the object's devi- 
ation from point-symmetry. While the atmospheric turbulence 
prevents us from measuring directly the phase, AMBER mea- 
sures two observables (DP and CP) which are related to the 
object's phase. DP approximately represents the object's phase 
in a spectral feature measured with respect to the continuurrfj 
Therefore, non-zero DPs mean a photocenter shift of the spec- 
tral feature forming region with respect to the continuum form- 
ing region. CP is the sum of phases around a closed triangle of 
baselines (i.e., <p\2 + f23 + fn) and not affected by the atmo- 
spheric turbulence. For point-symmetric objects, CP is always 
zero or n. Non-zero and non-7r CPs, whether in the continuum 
or in some spectral features, indicate asymmetry of the object. 
Two-telescope AMBER observations provide only one visibility 
and one DP (CP cannot be measured with two telescopes). 

Betelgeuse was observed on 2008 January 08 with AMBER 
using three ATs in the E0-G0-H0 linear array configuration with 
16-32^48 m baselines (AMBER Guaranteed Time Observation, 
Program ID: 080.D-0236A, PL: K. Ohnaka). Since the E0-G0- 
H0 configuration is a linear array lying at +7 1 from North (+90° 
= East) on the ground, the position angles of the three pro- 
jected baseline vectors are the same. We used the /f-band high- 
resolution mode (HR_K) of AMBER with a spectral resolution 
of 12000 covering wavelengths from 2.28 to 2.31 /mi. As shown 
below, this wavelength range was chosen to observe the strong 
12 C ie O (hereafter simply CO) first overtone lines near the (2,0) 
band head. The //-band brightness of Betelgeuse is too high for 
the VLTI fringe tracker FINITO. However, the extremely high 
brightness of Betelgeuse (K = -4.4), together with the excel- 
lent weather conditions (07 3-074 seeing), enabled us to detect 
low-contrast fringes on all three baselines without FINITO. 



1 Exactly speaking, two pieces of information are lost in the deriva- 
tion of DP from AMBER observations: the absolute phase offset and 
the linear phase gradient with respect to wavenumber. 
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We also downloaded AMBER data of Betelgeuse obtained 
on 2006 February 10 (Program ID: 60.A-9054A) from the ESO 
data archive. These data were taken with two ATs in the E0-G0- 
16 m configuration using the Zf-band medium-resolution mode 
(MR_K) without FINITO. The wavelength range between 2.1 
and 2.2 /mi was covered with a spectral resolution of 1500. 
Since there are no strong spectral features in this wavelength 
region, these MRJC data approximately sample the continuum. 
In both runs in 2008 and 2006, Sirius (a CMa, Al V, K = -1.4) 
was observed for the calibration of the interferometric data of 
Betelgeuse. We adopted an angular diameter of 5.6 mas for 
Sirius given by Richichi & Percheron (2005). Observations of 
stars with known angular diameters and no asymmetry are nec- 
essary to evaluate the so-called interferometer transfer function, 
which represents the instrumental and atmospheric effects on 
visibility and phase measurements, and to calibrate interferomet- 
ric data of a science target. A summary of the observations is 
given in Table Q] 

For the reduction of the AMBER data, we used amdlib 
ver.2.20, which is based on the P2VM algorithm (Tatulli et al. 
120071 1. We split each data set of Betelgeuse and Sirius into five 
or six subsets with each subset containing 500 frames and de- 
rived the visibilities, DPs, and CPs, as well as spectra, from each 
subset. One of the parameters in the reduction with amdlib is 
the frame selection criterion. For each subset, we checked for 
a systematic difference in the results by taking the best 20%, 
40%, 60%, 80%, and 100% of all frames in terms of the fringe 
S/N ratio. We found out that the visibilities obtained from the 
first two subsets of the Betelgeuse data #2 and the Sirius data 
taken in 2008 show a significant dependence on the selection 
criteria (the more frames are included, the lower the visibility 
becomes), while the visibilities from the other subsets do not 
show such a dependence. Furthermore, this dependence occurs 
only on the 16 m and 48 m (E0-G0 and E0-H0) baselines, while 
the visibilities on the 32 m (G0-H0) baseline are not affected by 
the selection criterion in any subset. The vibration of the AT at 
the E0 station is very likely to be the cause of this problem, be- 
cause it appears only on the baselines using the E0 station, and 
vibration always lowers visibility. Therefore, we dropped these 
subsets affected by the AT vibration for the derivation of the vis- 
ibilities, DPs, and CPs on the 16 m and 48 m baselines, while 
we used all subsets for the observables on the 32 m baseline. For 
the subsets not affected by the vibration, the reduction with dif- 
ferent selection criteria does not lead to a significant systematic 
difference in the results. The selection of fewer frames only re- 
sults in larger errors in the final results, while the inclusion of 
frames with very poor S/N ratios produces spurious results at 
a few wavelengths. Therefore, we included 80% of all frames in 
the subsets not affected by the vibration. For the AMBER MR_K 
data, we did not find a signature of vibration, so we used all sub- 
sets with 80% of frames included. Currently, AMBER data taken 
in the HRJC mode are affected by the Fabry-Perot effect caused 
by the InfraRed Image Sensor (IRIS), which stabilizes the image 
motion (see Fig. 6 in Weigelt et al. 2008). The IRIS Fabry-Perot 
effect is seen as time-dependent high-frequency beating in the 
raw visibilities, DPs, and CPs of Betelgeuse and Sirius plotted 
as a function of wavelength. However, fortunately, it is mostly 
removed by dividing the data of Betelgeuse with that of Sirius, 
and the IRIS beating is barely discernible in the calibrated visi- 
bilities, DPs, and CPs of Betelgeuse. This is because the data of 
Betelgeuse and Sirius were taken close in time. In other cases, 
the IRIS beating degrades the final data significantly. 



2 Available at http://www.jmmc.fr/data_processing_amber.htm 



While the visibilities and DPs on the shortest baseline re- 
duced from the data taken with a spectral resolution of 12000 are 
of sufficient quality, the visibilities and DPs on the longer base- 
lines and CPs turned out to be noisy. Therefore, we improved the 
S/N ratio for these observables by binning the data (object, dark, 
sky, and P2VM calibration data) in the spectral direction. For the 
visibilities and DPs on the middle baseline, binning with a box 
car filter with a width of three pixels turned out to be sufficient. 
This results in a spectral resolution of 8000 instead of 12000 
achieved with the two-pixel sampling. For the observables on 
the longest baseline and CPs, it was necessary to bin the data 
with five pixels, corresponding to a resolution of 4800, to obtain 
reasonable S/N ratios. As shown below, the individual CO lines 
can still be resolved with these lowered spectral resolutions. 

The errors of the resulting visibilities, DPs, and CPs were es- 
timated from the standard deviation among the results obtained 
from five or six subsets. The errors of the calibrated observables 
were computed from such errors in the data of Betelgeuse and 
Sirius. Since we have only one data set of Sirius, it is impos- 
sible to estimate the systematic error in the transfer function. 
Therefore, to account for this error source, we added a system- 
atic error of 5% to the above errors. 

In both runs in 2008 and 2006, Sirius served not only as an 
interferometric calibrator but also as a spectroscopic standard 
star. We attempted to remove telluric lines as much as possi- 
ble by dividing the spectra of Betelgeuse with that of Sirius, 
although the difference in air mass for the 2008 data did not 
allow us to achieve this perfectly. The telluric lines identified 
in the spectrum of Sirius were used for wavelength calibration. 
As a template of the telluric lines, we convolved the atmospheric 
transmission spectra from Wallace & Hinkle d 19961 1 to match the 
resolutions of the HRJC and MR_K modes of AMBER. The un- 
certainty in wavelength calibration is ~6 X 10~ 5 /mi (~8 km s _1 ) 
for the HRJC observations and ~5xl0~ 4 /mi (~70km s _1 ) for the 
MRJC observations. The wavelength scale was then converted to 
the heliocentric frame using the IRAr0 task RVCORRECT. 

2.2. VINCI 

To discuss temporal variations in the Zf-band visibility, we 
also downloaded interferometric data of Betelgeuse taken with 
VINCI (VLT INterferometer Commissioning Instrument) from 
the ESO archive (Program ID: 60.A-9222A). These data were 
obtained as part of the commissioning of the instrument. A de- 
tailed description of the instrument is given in Kervella et al. 
(2000). As summarized in Table 12 the VINCI observations of 
Betelgeuse occurred on eight nights between 2001 and 2003 us- 
ing two 40 cm siderostats in the E0-G0-16m and B3-C3-8m con- 
figurations (both configurations lie at +71° from North (+90° = 
East) on the ground). A number of calibrators were observed on 
these nights (6-35 calibrator measurements on each night), as 
listed in Table[3] 

We used the VINCI data reduction software ver.3.C0 
(Kervella et al. 120041 ) to derive visibility. The interferometer 
transfer function was computed from all calibrator measure- 
ments taken during a given night, and the mean of these trans- 
fer function values was used to obtain the calibrated visibili- 
ties of the science target. The error of each calibrated visibil- 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

4 Available at http://www.jmmc.fr/data_processing_vinci.htm 



4 



Ohnaka et al.: Spatially resolving the inhomogeneous atmosphere of Betelgeuse 



ity was derived from the statistical error in each measurement 
of the science target given by the reduction software and the 
error in the transfer function. This latter error results from the 
statistical error in each calibrator measurement and the standard 
deviation of the transfer function values obtained on the given 
night. The VINCI data reduction software computes visibility 
using two different algorithms: Fourier transform and wavelet 
transform. For the data presented here, the calibrated visibil- 
ities derived with two methods agree well, and we only give 
the results obtained with the wavelet transform in Table [2] The 
errors in the calibrated visibilities are typically 1-3%. Unlike 
AMBER, VINCI observations were made with the /^-broadband 
filter covering from 2 to 2.4 /mi. Therefore, we computed the 
effective wavelength using the VINCI transmission presented 
in Wittkowski et al. (2004) and the spectrum of Betelgeuse 
observed with the Stratoscope II (detector B) by Woolf et al. 
( 119641 1. The resulting effective wavelength, 2.175 /mi, was used 
for calculating the spatial frequency for each observation. 



3. Results 

FigureQ]shows the calibrated visibilities, DPs, and CPs observed 
toward Betelgeuse as a function of wavelength. The visibilities 
and DPs on the middle and longest baselines, as well as the CPs 
were derived from the binned data (spectral resolutions of 8000 
and 4800), while the results on the shortest baseline were de- 
rived from the data with a spectral resolution of 12000. The fig- 
ure reveals the detection of interferometric signals even on the 
longest baseline (48 m). This marks the highest spatial resolution 
(9 mas) obtained for Betelgeuse, corresponding to nearly five 
resolution elements over its stellar disk. The visibilities, DPs, 
and CPs derived from the two data sets show mostly the same 
spectral features, demonstrating that the observed features are 
real despite the low fringe contrast. There seems to be a differ- 
ence in visibility level for the E0-G0-16m baseline. However, 
given the uncertainties in the calibrated visibilities shown in 
the figure, this discrepancy is marginal. We note that the CPs 
near the CO band head between 2.294 and 2.296 /mi are much 
noisier — errors of 50-100° — than at the other wavelengths even 
in the binned data. 

3.1. Continuum: A < 2.293 /<m 

The observed spectrum of Betelgeuse below 2.293 /mi shows 
only several weak atomic and molecular absorption features, as 
identified in Fig. There are also subtle signatures in the ob- 
served visibilities corresponding to these features. In particular, 
the effects of the Ti (+ HF) feature at 2.29 //m can be seen in the 
visibilities on all three baselines and also possibly in the CPs. 
These weak lines form in the deep photospheric layers and can 
be used for testing photospheric models, which we plan in a fu- 
ture paper. 

The spectral resolution of the HR mode of AMBER allows 
us to select continuum points which are not affected by these 
lines. Since the visibilities on the longest baseline were derived 
with the five-pixel binning, we also derived the visibilities on 
the shorter baselines and the spectra from the data binned with 
five pixels. Then we selected 37 continuum points in the spectra, 
avoiding the lines. In Fig. [2] we plot the visibilities observed at 
these continuum points below 2.293 /mi as a function of spa- 
tial frequency. Also plotted are the AMBER MR_K data be- 
tween 2.1 and 2.2 /mi, VINCI data, /f-band data obtained at the 
Infrared Optical Telescope Array (IOTA) by Perrin et al. (2004), 



and the /f-band measurements by Dyck et al. (1992). The vis- 
ibilities derived from the AMBER MR_K data agree very well 
with the result obtained by Hernandez & Chelli (120071 ) using 
not only amdlib but also another algorithm based on the Fourier 
transform. The errors in the visibilities from the MR_K data are 
~6%. The DP between 2.1 and 2.2 /mi derived from these data 
is zero within a measurement error of ~5°, which is no surprise, 
given the absence of strong spectral features in this wavelength 
range. While the AMBER HR_K and MR_K data measure the 
continuum almost free from the effects of molecular/atomic fea- 
tures, the VINCI and IOTA data as well as the measurements of 
Dyck et al. ( 119921 ) were taken with a broad band filter spanning 
the entire K band, which includes the molecular absorption fea- 
tures due to CO, H2O, and CN and atomic lines. However, since 
the CO and H2O bands appear only at limited wavelengths (CO: 
longward of 2.3 /mi, H2O: either edge of the K band), and CN 
and atomic features are weak in Betelgeuse, the total /f-band 
flux is dominated by the continuum. Therefore, the data taken 
with the K broad-band filter represent the visibilities in the con- 
tinuum in first approximation, and this is why they are included 
in the discussion below. 

To derive the angular size of the object in the continuum, we 
fitted the observed visibilities with a uniform disk and a limb- 
darkened disk. Uniform-disk fitting to the AMBER continuum 
data in 2008 and 2006 results in 43.19 + 0.03 mas (reduced x 2 
= 4.6) and 42.69 + 0.01 mas (reduced^ 2 = 0.01), respectively, 
while the fitting to all the data (AMBER, IOTA, and VINCI) 
results in 43.16 + 0.03 mas (reduced^ 2 = 3.3). Fitting to the 
AMBER HR_K data with a power-law-type limb-darkened disk 
(Hestroffer et al. 119971 ) results in a limb-darkened disk diam- 
eter of 43.56 + 0.06 mas and a limb-darkening parameter of 
(1.2 + 0.07) x 10" 1 (reduced x 2 = 3.44), while fitting to all the 
data results in a limb-darkened disk diameter of 43 .50 + 0.04 mas 
and a limb-darkening parameter of (9.0 + 0.4) x 10~ 2 (reduced 
X 1 = 2.5). These results agree well with the limb-darkened disk 
diameter and the limb-darkening parameter derived by Perrin et 
al. (f2004) from the ^-broadband IOTA data. Since a uniform 
disk has a limb-darkening parameter of 0, the fit suggests a small 
limb-darkening effect in the ^-band continuum. Limb-darkened 
disk fitting to the AMBER MR_K data did not give a meaningful 
result, because the data are located in one narrow spatial fre- 
quency range. The fit to the AMBER HR_K data in 2008 shows 
that the 16, 32, and 48 m baselines correspond to the second, 
third, and fifth lobes. Although the longest baseline falls near 
the fourth null (i.e., between the fourth and fifth lobes), the de- 
tection of fringes in the fifth lobe is corroborated by the 180° 
CPs measured in the continuum below 2.293 /mi (Fig. QJ). This 
is because CP is the sum of the phases on the three baselines 
(E0-G0, G0-H0, and H0-E0), that is, CP = 180° (2nd lobe) + 0° 
(3rd lobe) + 0° (5th lobe) = 180°. While the above fit shows the 
AMBER HR_K data cannot be perfectly fitted with a uniform 
disk or limb-darkened disk (i.e., x 2 > IX Fig- E] suggests that 
the deviation from the fitted curves is not drastic at the time of 
the observations. This means that the spatial scale of inhomo- 
geneities in the A'-band continuum is even smaller than the res- 
olution of the longest baseline (9 mas) and/or the contrast of the 
inhomogeneities is small. This picture is also supported by the 
result that the CPs measured in the continuum are 180° within 
the errors, as expected from the above fit with a uniform disk or 
a limb-darkened disk. 

The data in the second and third lobes were taken at seven 
different epochs as given in Fig. [2] Therefore, these data points 
contain information about the effects of time-dependent sur- 
face inhomogeneities on the ^T-band visibility in the contin- 
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uum. Recent 3-D hydrodynamical convection simulations for 
RSGs by Chiavassa et al. ( 120071 and Ludwig & Beckers ( f2008T > 
predict significant temporal variations in visibility beyond the 
first null due to ever-changing inhomogeneities. For example, 
Chiavassa et al. (12007ft present the amplitude of such variations 
in the 2.2 fim visibility in the second and third lobes predicted 
for Betelgeuse, which are plotted by the dotted lines in Fig. [2] 
(the model visibility at 2.1 //m presented in Ludwig & Beckers 
(2008 ) also shows a similar behavior). However, most of the ob- 
servational data points in the second and third lobes are located 
near the upper boundary or at least in the upper half of the pre- 
dicted range, instead of equally distributed above and below the 
center of the range. This means that the model predicts the vis- 
ibility to be systematically lower than the observations. The ob- 
served data in the second lobe, where we have AMBER, VINCI, 
and IOTA data from five epochs, show only a modest scatter 
around the limb-darkened-disk fit. The data on the third lobe do 
not show large variations, either, although we have data from 
only two epochs. Therefore, the AMBER, VINCI, and IOTA ob- 
servations imply that the 3-D convection simulation of Chiavassa 
et al. (120071) predicts too large deviations (presumably due to too 
strong inhomogeneities) from the limb-darkening disk which fits 
the observed data reasonably in the /f-band continuum. 

However, given that only our AMBER HR_K data probes the 
very high-spatial frequency regime and the number of the obser- 
vation epochs is not yet very large, we cannot conclude whether 
or not the visibility at these high spatial frequencies shows lit- 
tle deviation from a uniform-dirk or limb-darkened disk all the 
time. Also, only our AMBER observations have measured CP, 
whose deviations from zero or n would represent a clear signa- 
ture of asymmetry. The visibilities of Betelgeuse observed with 
the //-broadband filter with IOTA show more deviation from a 
uniform disk than in the K band continuum, and non-zero CPs 
are also detected in the third and fourth lobes (Haubois et al. 
2006). While these results reveal inhomogeneous structures seen 
with the //-broadband filter (including the absorption features 
due to molecules such as H2O, CO, and OH), it is not yet clear 
whether the amplitude of the deviation is consistent with the cur- 
rent convection simulation models, because the //-band IOTA 
data represent only one epoch. Further spectro-interferometric 
monitoring observations in the K and H bands at such high spa- 
tial frequencies as obtained here will enable us to study the spa- 
tial and time scale of inhomogeneities in the continuum and pro- 
vide stronger observational tests for the 3-D simulations of con- 
vection in RSGs. Particularly, multi-epoch observations at the 
same spatial frequencies will provide a direct test for the visibil- 
ity fluctuation predicted by the simulations. 

Recently, Townes et al. (2009) have reported a monotonic 
decrease by 15% in the 1 1 fim diameter of Betelgeuse in the past 
15 years, from 1993 to 2009. Its origin is unclear. We examined 
a possible long-term variation in the Zf-band uniform-disk diam- 
eter in the interferometric observations in the literature and the 
VINCI data. The uniform-disk diameters derived from the first- 
lobe data in the past are 44.2 + 0.2 mas (Aug. -Sep. 1990, Dyck 
et al. [19921 ) and 43.26 + 0.04 mas (Nov. 1996, Perrin et al. |2004l 
We fitted the VINCI data in the first lobe taken in Dec. 2002 and 
Jan. 2003 with a uniform disk model computed with the Zf-band 
transmission of VINCI and the spectrum of Betelgeuse, which 
were used for computing the effective wavelength. The resulting 
uniform-disk diameter, 42.02 ± 0.4 mas, agrees very well with 
the 42.06+0.4 obtained by Meisner (priv. comm.) from the same 
VINCI data but using a different algorithm based on coherent 
integration (Meisner 120031 ). The error includes the uncertainty 
in the determination of the effective wavelength in VINCI ob- 



servations (Meisner, priv. comm.), which we also added to the 
error in our uniform-disk diameter. Therefore, the Zf -broadband 
uniform-disk diameter shows a decrease by 5 + 2% from 1990 
to 2003, much less pronounced than at 1 1 fim. The continuum 
uniform-disk diameters of 42.69 and 43.19 mas obtained from 
the AMBER data taken in 2006 and 2008 might indicate an in- 
crease in the angular size after 2003. However, these diameters 
obtained from the data beyond the first null can be affected by 
inhomogeneities (although not very strong), which makes the 
apparent increase in the angular size inconclusive. 

3.2. CO first overtone lines: A > 2.293 fim 

Figure Q]reveals salient signatures of the CO first overtone lines 
in the visibilities, DPs, and CPs longward of 2.293 fan. Figure|3] 
shows the enlarged views between 2.299 and 2.306 fim (only the 
data set #1 is plotted for visual clarity), where the observed spec- 
trum and the line positions are overplotted in each panel to show 
the shape of the visibilities, DPs, and CPs within the CO lines. 
Each absorption feature consists of two transitions with high and 
low rotational quantum numbers /. But for simplicity, we refer 
to such an absorption feature with high and low J as a "line" in 
the present work. Comparison of the observed line positions with 
the laboratory data results in a heliocentric velocity of 24 km s _1 . 
This value agrees with the velocity range 17-27 km s _1 derived 
from the visible and IR atomic and molecular lines as well as the 
m m CO lines (e.g., Bro oke et al. 119741 Huggin s fT987l Smith et 
al. [19891 Huggins et al. [T994l Ryde et al. llfj?9l . 

Figure [3^ shows the visibilities observed on the three base- 
lines. Particularly surprising is that the visibility observed within 
a given CO line on the shortest baseline (~16 m) is anti- 
symmetric with respect to the line core (i.e., "-"-shaped). While 
the visibility on the middle baseline (~32 m) is roughly symmet- 
ric, the visibility on the longest baseline (-48 m) is asymmetric 
with the peak slightly redshifted with respect to the line core in 
most cases. These results mean that Betelgeuse appears different 
in the blue and red wings of the CO lines and that the blue and 
red wings originate in spatially distinct regions differing in size 
and/or shape. The observed DPs and CPs show remarkable non- 
zero and non-7r values, as large as DP = -130° (Fig.[Tg) or CP = 
180°+ 90° (Figs.QJ and[3j;). These non-zero and non-7rDPs and 
CPs confirm that the blue and red wings of the CO lines originate 
in spatially distinct regions. On the other hand, the observed DPs 
and CPs are nearly zero near the line core, which means that the 
star appears symmetric. 

One might suspect that such blue-red asymmetry in visibili- 
ties, DPs, and CPs within the CO line profiles can be explained 
by the fact that one CO absorption feature, which appears to be 
a single line, is a blend of two transitions with low and high /. If 
two lines with different excitation potentials form at different re- 
gions over the stellar surface, it may make the star appear differ- 
ent in the blue and red wings and cause the blue-red asymmetry. 
However, the relative positions of the high and low J lines swap 
at 2.3032 fim: shortward of this wavelength, the high J lines ap- 
pear blueward of the low J lines, while above 2.3032 fim, the 
high J lines appear redward of the low J lines. Still, the same 
asymmetry is observed in all CO lines, whether shortward or 
longward of 2.3032 fim. Therefore, the observed asymmetry in 
the visibility, DPs, and CP cannot be explained by the blend of 
two transitions with high and low J. 

Stellar rotation or spherically expanding/infalling flows are 
also unlikely to be the cause of the blue-red asymmetry in visi- 
bilities and phases. The projected photospheric rotational veloc- 
ity of Betelgeuse is small, vsin; w 2-5 kms' 1 (Uitenbroeck 
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1998; Harper & Brown 120061 ) . Moreover, the rotational axis 
inferred from the chromospheric emission lines lies at ~65° 
(Harper & Brown 2006 ), which is very close to the position angle 
of the projected baselines of our AMBER observations. In this 
case, rotation can have no noticeable effects on visibilities and 
phases. In spherically expanding/infalling flows, the velocity in 
the line of sight changes as a function of the angular distance 
from the stellar disk center. Therefore, the photons at different 
wavelengths (= different velocities) within a line profile origi- 
nate in annular regions with different sizes, which makes the star 
appear different in the blue and red wings. However, obviously, 
such spherically symmetric models cannot explain the observed 
non-zero/non-7r DPs and CPs. 

The observed blue-red asymmetry in visibility, DP, and CP 
may be explained by an inhomogeneous velocity field in the at- 
mosphere, in which upwelling and downdrafting CO gas exists 
in spatially distinct regions. In spectroscopic analyses using 1-D 
model atmospheres, such a non-thermal velocity field is empiri- 
cally incorporated as macroturbulence, which manifests itself as 
a broadening of spectral lines ("macro" means the spatial scale 
of the non-thermal gas motion is larger than the length of the 
unit optical depth of photons). In the next section, we examine 
whether the observed visibilities, DPs, and CPs can be explained 
by such an inhomogeneous velocity field. 

4. Modeling 

To characterize the inhomogeneous velocity field and the prop- 
erties of the CO gas in the atmosphere of Betelgeuse, we con- 
structed the following patchy two-layer model, in which the star 
is surrounded by the inner and outer CO layers. The star was as- 
sumed to be a blackbody of an effective temperature of 3600 K 
based on the value derived by Perrin et al. ( 120041 1. The inner 
layer represents the photosphere — the region included in 1-D 
photospheric models, usually spanning from continuum optical 
depths of ~ 1 1 — 1 2 to ~ 10~ 5 -10~ 6 . The outer layer represents 
the MOLsphere. Given the presence of dense H2O gas in the 
MOLsphere suggested for Betelgeuse from spectroscopic and 
interferometric observations (Perrin et al. 2004, 2007; Ohnaka 
2004a] Tsuji 120061 ), it is plausible that there is also a significant 
amount of CO in the MOLsphere. 

The basic picture of our model is depicted in Fig. [4] This 
two-layer model is similar to that described in Ohnaka (2004b), 
but we introduced two modifications. Firstly, the geometrical 
thickness of each layer is assumed to be very small compared to 
its radius, and therefore, the two layers are detached from each 
other. This simplifies the computation of the line opacity in the 
presence of a velocity field as described below. Secondly, while 
the column density and temperature of each layer were assumed 
to be constant over the stellar surface, we incorporated the fol- 
lowing inhomogeneous velocity field: CO gas is moving radially 
outward (or inward) with a velocity of Vfl OW m one patch, while 
it is moving in the opposite direction at the same velocity in the 
remaining region. We assumed the same velocity field for two 
layers. As Fig. |4j? illustrates, we define such a patch as a cone. 
Its half-opening angle (0) characterizes the size of the patch, 
while its position is specified by two angles, 9 and <f>, which de- 
fine the vector connecting the center of the patch and the center 
of the star. 

We do not know a priori the actual number and shape of 
patches (if any) on the stellar surface, and the present data are in- 
sufficient to derive the inhomogeneous surface pattern uniquely. 
We assume only one patch in our modeling to keep the number 
of free parameters as small as possible. The aim of our modeling 



is not to derive the actual inhomogeneous surface structure but 
to see if there is indeed a model with an inhomogeneous veloc- 
ity field which can explain the observed blue-red asymmetry in 
visibility, DP, and CP within the CO lines. 

To decrease the number of free parameters, we fixed the 
CO column density and radius of the inner layer as follows. A 
spherical photospheric model with the molecular opacities due to 
CO, TiO, H2O, OH, and SiO incorporated based on the opacity 
sampling (Ohnaka, in prep.) was computed using the stellar pa- 
rameters of Betelgeuse given in Tsuji (2006). This photospheric 
model gives a CO column density of 5 x 10 22 cm 2 , which we 
adopted for the inner layer. We assumed the radius of the inner 
CO layer to be a half of the geometrical thickness of this pho- 
tospheric model (~0.1 R*), that is, a radius of 1.05 R*. While 
this choice is ambiguous, it turned out not to affect the result 
significantly, as far as the radius is smaller than ~1.1 It is 
necessary to include the microturbulence, which represents the 
non-thermal gas motion on a spatial scale smaller than the length 
of the unit optical depth of photons. The microturbulent veloc- 
ity in the photosphere of Betelgeuse derived from spectroscopic 
analyses ranges from 4 to 6 km s -1 (Lambert et al. 1 19841 Tsuji et 
al. 119941 Tsuji 2006). We adopted a microturbulent velocity of 
5 km s _I for both the inner and outer CO layers in our modeling. 

For the outer layer, we fixed its radius and temperature based 
on the previous spectroscopic and interferometric studies of 
the MOLsphere. The radius of the MOLsphere of Betelgeuse 
measured from the near- and mid-IR H2O features is 1.45 R* 
(Ohnaka 2004a), 1.3 R+ (Tsuji |2006j ), and 1.31-1.43 fl* (Perrin 
et al. 120071 ) . We tentatively adopted 1.45 R+ for our model- 
ing. The temperature of the H2O MOLsphere derived by the 
above authors range from 1500 to 2250 K. A gas temperature 
of 1800 K, which lies roughly in the middle of this range, was 
adopted for the outer CO layer. Therefore, the free parameters 
of our model are the temperature of the inner CO layer (Tin), the 
CO column density of the outer CO layer {N out ), the velocity of 
the CO gas motion (vg ow ), and the position and size of the patch 
(6, 4>, and 0). 

The intensity distribution for this patchy model was com- 
puted at each wavelength between 2.291 and 2.309 //m by per- 
forming ray tracing along a number of lines of sight (see Ohnaka 
2004a 1 2004b ). The only difference from these previous studies 
is the inclusion of the Doppler shift due to the velocity field in 
the calculation of the line opacity. Since each layer is geometri- 
cally thin, the velocity component along the line of sight within 
one layer can be approximated to be constant, which simplifies 
the computation. The CO line list of Goorvitch ( 1994 ) was used 
for the calculation of the line opacity. The monochromatic 2-D 
complex visibility was calculated by taking the Fourier trans- 
form of the intensity and then spectrally convolved to match 
the resolution of our AMBER observations (12000, 8000, and 
4800). The visibility amplitude and phase were derived from this 
spectrally convolved 2-D complex visibility, which was scaled 
with the uniform-disk diameter of 43.19 mas derived from the 
AMBER HR_K data in the continuum. CP was obtained as the 
sum of Fourier phases predicted for three baselines, while DP 
on a given baseline was computed by taking the difference in 
Fourier phase between in a CO line and in the continuum below 
2.293 yum. The correspondence between the sign of CP measured 
with AMBER and the positional offset in the plane of the sky is 
described in Kraus et al. (2009]). 

Figure [5] shows a model with T m = 2250 K and N out 
I x 10 20 cirT 2 , characterized by a large, upwelling spot (0 = 
60°, Vfl ow = 10 kms -1 ) covering nearly a half of the apparent 
stellar disk (6 = 40°, 4> = 10 °, see also Figs.|6h and|6}l). Given 
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the simple nature of the model, the overall agreement is reason- 
able, although there are still discrepancies between the model 
and the observed data. Particularly, the agreement is poor near 
the band head, which is discussed below. Also, the visibility on 
the longest baseline as well as the CP predicted by the model is 
too smooth compared to the observed data (note that the spec- 
tral resolution of the model matches that of the binned data). 
Figure [6] shows the images predicted by this model at different 
wavelengths within a CO line (R(29) and R(72)). The figure il- 
lustrates that the velocity field makes the star appear different 
in the blue and red wings. This can cause the visibility to be 
anti-symmetric with respect to the line core as observed on the 
shortest baseline and also explains the observed asymmetric DPs 
and CPs. The stellar image at the line core is the sum of the blue- 
and red-shifted components. This makes the star appear nearly 
centrosymmetric (Fig. [6}:) and results in DPs and CPs close to 
zero at the line core as observed. 

The agreement between the observed data and the model be- 
comes poorer near the band head. In particular, the deviations 
of DP from zero predicted for the middle and longest baselines 
are too small compared to the observed values. The reason why 
the DPs predicted by the model do not show significant devia- 
tion from zero is as follows. Since the lines are very crowded 
near the band head, the blue (or red) wing of one line overlaps 
with the red (or blue) wing of the adjacent lines. This means 
that the appearance of the star at a given wavelength in such a 
crowded region is the sum of the blue- and red-shifted compo- 
nents as in the case of the image at the core of an isolated line. As 
Fig. [6J; shows, the resulting image appears roughly centrosym- 
metric. This is the reason why the DPs predicted by the model 
do not show noticeable deviations from zero. A possible reason 
for the discrepancy near the band head is our assumption that 
the CO column density and temperature are uniform over the 
stellar surface. Inhomogeneities in CO column density and/or in 
temperature over the surface may reconcile the disagreement be- 
tween the model and the observed data near the band head. 

The uncertainty ranges of T\ n , N out , and vo ow derived from 
our modeling are 2000-2500 K, 5 x 10 19 -2 x 10 20 cnr 2 , and 
10-15 km s , respectively. We also computed models with 7V out 
set to zero (i.e., no MOLsphere), but such models cannot repro- 
duce the visibilities observed near the band head. The uncer- 
tainties in the position and size of the patch are very large: the 
ranges for 0, 6, and (p are 20°-80°, 20°-80°, and 0°^0°, re- 
spectively. However, these uncertainties in the position and size 
of the patch should not be taken at face value, because we as- 
sumed for simplicity that there is only one spot, which may not 
be the case. Also, we cannot conclude whether the CO gas in 
the patch is moving upward or downward. We computed mod- 
els in which the patch is moving downward and found out that 
there are parameter sets which are consistent with the observa- 
tions. However, apart from these uncertainties in the actual sur- 
face pattern and the direction of the velocity field, our AMBER 
data and modeling suggest an inhomogeneous velocity field with 
amplitudes of ~ 1 0— 1 5 km s _1 . 

5. Discussion 

The amplitude of the velocity field suggested from our mod- 
eling, 10-15 km s , compares favorably with the macroturbu- 
lent velocities derived in the previous spectroscopic analyses. 
Macroturbulence manifests itself as the broadening of spectral 
lines in observed spectra. In spectral analyses using 1-D model 
atmospheres, macroturbulence is usually incorporated as an ad- 
ditional line broadening factor with a Gaussian distribution (i.e. 



°c exp(-(/l//l macro ) 2 )). The macroturbulent velocity — either as 
the Gaussian dispersion V macro or the FWHM = 1.665V macro cor- 
responding to A m - dcm — is derived so that synthetic spectra repro- 
duce observed line profiles. For Betelgeuse, macroturbulent ve- 
locities as high as 20 km s -1 (FWHM) were derived from the 
optical lines by Lobel & Dupree (2000) and Gray (120001) . while 
smaller values of 10-12 km s _1 were obtained from near- and 
mid-IR lines (Jennings et al. 1986; Jennings & Sada 1998; Ryde 
et al.|2006a| Tsuji|2UU6]>. 

Josselin & Plez (2007) analyzed spectral line profiles in the 
optical and extracted information about the velocity field in a 
sample of RSGs. For Betelgeuse, they found that the strong lines 
with lower excitation potentials of ~1 eV show blue- and red- 
shifted components at approximately ±10 kms~' with respect 
to the velocity of weak lines with excitation potentials of ~3 eV. 
These velocities roughly agree with the 10-15 km s _1 derived 
from our modeling of the CO first overtone lines with excitation 
potentials of <1 eV. Josselin & Plez (2007) detected no peri- 
odic or regular temporal variation in the velocity of the blue- 
and red-shifted components, which led them to interpret these 
two components as upward and downward convective motion. 

On the other hand, it is not so obvious whether or not such 
strong convective motion can be present in the CO line formation 
layers, where the convective energy flux is small. Alternatively, it 
is possible that the CO gas motion detected by our AMBER ob- 
servations may represent local mass ejections. High-resolution 
studies of the dusty RSG VY CMa and less massive evolved 
stars such as IRC+10216 and CIT6 suggest that the mass loss 
is accompanied by episodic mass ejections in clumps or arcs 
(Humphreys e t al. [20071 Smith et al |2009l Weigelt et al. l2002l 
Monnier et al. 2000). While the circumstellar envelope around 
Betelgeuse is spherical on a large scale unlike VY CMa, inho- 
mogeneities such as clumps and plumes have been detected (Plez 
& Lambert 120021 Smith et al. |2TX)9"l ). The long-term variability 
of the Ha line velocity can also be interpreted as a consequence 
of "intermittent failed ejections", in which material is flung out 
and falls back toward the star (Smith et al. |1989l ). 

To summarize, our AMBER observations are the first 
spatially-resolved detection of macroturbulent gas motion in a 
stellar atmosphere (photosphere and possibly MOLsphere) other 
than the Suifl The spatially resolved CO gas motion can repre- 
sent the convective motion in the photosphere (and also in the 
MOLsphere) or the motion related to intermittent mass ejections 
in clumps and/or arcs. 

The (u, v) coverage of the present AMBER data is insuffi- 
cient to derive the actual inhomogeneous structure in the pho- 
tosphere and MOLsphere. In particular, the lack of visibilities 
in the first lobe makes it difficult to measure the size of the 
MOLsphere seen in the CO lines, and we simply assumed the 
same radius and temperature as those derived from the near- and 
mid-IR H2O features. While baselines shorter than the 16 m cur- 
rently available at the VLTI is desirable, it is also possible to 
observe Betelgeuse in the first lobe by taking advantage of the 
projection effect. Such new data will be indispensable for con- 
straining the geometrical extent of the CO gas in the MOLsphere 
of Betelgeuse. 

Lastly, we estimate the number of observations necessary to 
reconstruct an image. As a rule of thumb, the number of (w, v) 
points should be larger than the number of filled pixels (i.e., 
pixels with stellar flux) in the reconstructed image (e.g., Haniff 



5 The chromospheric gas motion toward Betelgeuse was spatially re- 
solved with the HST observations in the UV by Gilliland & Dupree 
fT996l >. 
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120071 1. This means that we need at least five (u, v) points — 
realistically about 10 points — to reconstruct a 1-D image with 
five resolution elements over the stellar disk as in our present 
work. Therefore, to reconstruct a 2-D image with 5x5 pixels, we 
need approximately 100 (u, v) points, which corresponds to ~33 
observations (~3 nights) with AMBER using three telescopes. 
It is crucial that these (u, v) points are as uniformly distributed 
as possible, which will become feasible when more new VLTI 
configurations, particularly short baselines perpendicular to E0- 
G0-H0, are opened. The above estimate is roughly consistent 
with the image reconstruction simulation for the VLTI 2nd gen- 
eration instrument MATISSE (Hofmann et al. 120081 . 

6. Concluding remarks 

We have spatially resolved the CO gas motion in the atmosphere 
of Betelgeuse with high-spectral resolution using VLTI/AMBER 
and successfully measured the visibilities, DPs, and CPs in the 
second, third, and fifth lobes, marking the highest spatial resolu- 
tion (~9 mas) achieved for Betelgeuse. The visibilities observed 
in the CO first overtone lines suggest that the blue and red wings 
of individual lines originate in spatially distinct regions, and the 
non-zero/non-7r DPs and CPs observed in the CO lines corrobo- 
rate this picture. Our simple model suggests an inhomogeneous 
velocity field with amplitudes of ~10-15 kms 1 in the atmo- 
sphere of Betelgeuse. These AMBER observations are the first 
to spatially resolve the so-called macroturbulence in a stellar at- 
mosphere other than the Sun. The spatially resolved CO gas mo- 
tion is likely to correspond to the convective motion in the upper 
photosphere (and possibly MOLsphere as well) or intermittent, 
failed clumpy mass ejections. Our modeling also shows that the 
AMBER data are consistent with the presence of the MOLsphere 
extending to ~ 1 .45 R+ with a CO column density of ~ 10 20 cm -2 
and a temperature of 1800 K. 

The visibilities and CPs observed in the continuum below 
2.293 yum do not show a drastic deviation from a limb-darkened 
disk with a diameter of 43.56 mas and a limb-darkening pa- 
rameter of 0.12 at the time of the observations. Comparison of 
our AMBER HR_K data with the previous AMBER MR_K data, 
VINCI data, and IOTA measurements reveals that the recent 3-D 
convection simulations for Betelgeuse predict the Tif-band con- 
tinuum visibility to be too low beyond the first null. 

We plan to continue AMBER observations of Betelgeuse to 
study temporal variations in the visibilities and phases in the CO 
lines as well as in the continuum. Such high-spatial resolution 
data will provide tight constraints on the time scale and spatial 
scale of the inhomogeneities. When more baselines, particularly 
short ones, become available, reconstruction of a 5 x 5-pixel 
image will be feasible with observations in ~3 nights. High- 
resolution imaging for various molecular lines will be essential 
for understanding the physical processes responsible for the in- 
homogeneous structures in RSGs. 
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Fig. 1. AMBER observations of Betelgeuse. In all panels except for a, the black and red lines represent the data set #1 and #2, 
respectively. The spectrum shown in the panel a was derived from the merged data. The positions of the CO first overtone lines as 
well as other atomic and molecular lines are marked with the ticks. In the panels except for a, two error bars near the left ordinate 
represent the typical errors in the continuum (left) and in the CO lines (right) for the data set #1. The errors for the data set #2 are 
shown near the right ordinate in the same manner. The error in the normalized spectrum is 0.5% and 1 % in the continuum and in the 
CO lines, respectively. The wavelength scale is in the heliocentric frame, a: Normalized flux, b-d: Visibilities observed on the E0- 
G0-16m, G0-H0-32m, and E0-H0-48m baselines. The visibilities on the middle and longest baselines are binned with a box car filter 
with widths of three and five pixels, respectively, e: Closure phases with five-pixel binning. Note that the errors in CP near the CO 
band head between 2.294 and 2.296 //m are as large as 50-100°. f-h: Differential phases observed on the E0-G0-16m, G0-H0-32m, 
and E0-H0-48m baselines. The DPs on the middle and longest baselines are binned with three and five pixels, respectively. 
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Fig. 2. /f-band continuum/broadband visibilities of Betelgeuse plotted as a function of spatial frequency. The insets show enlarged 
views of the second, third, and fifth lobes. The error bars of the single AMBER data points are exemplarily shown in the insets. The 
errors of the VINCI and IOTA data are also shown in the insets. The solid and dashed lines represent the visibilities for a uniform 
disk with a diameter of 43.19 mas and for a limb-darkened disk with a diameter of 43.56 mas and a limb-darkening parameter of 
0.12 (power-law-type limb-darkened disk of Hestroffer (1997 ), respectively. The dotted lines represent the full amplitude of the 
variations in the 2.22 /jm visibility due to time-dependent inhomogeneous structures predicted by the 3-D convection simulation of 
Chiavassa et al. d20071 l. 
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Fig. 3. Visibilities (a), differential phases (b), and closure phases 
(c) observed in the CO first overtone lines toward Betelgeuse. 
In the panels a and b, the black, red, and blue lines represent 
the visibilities or differential phases observed on the E0-G0- 
16m, G0-H0-32m, and E0-H0-48m baselines, respectively. In 
the panel a, the visibility obtained on the longest baseline (blue) 
is scaled by a factor of six and shifted downward by 0. 1 for the 
sake of visual clarity. The black, red, and blue dotted lines repre- 
sent the continuum visibilities of a uniform disk with 43. 19 mas 
for the E0-G0-16m, G0-H0-32m, and E0-H0-48m baselines, re- 
spectively. In the panel c, the observed CP is shown by the black 
solid line. The dotted line represents CP = n, which is observed 
in the continuum. The visibilities and DPs on the middle and 
longest baselines were derived from the data binned with three 
and five pixels, respectively. The CP was derived from the data 
binned with five pixels. In each panel, the normalized flux (with- 
out binning) is overplotted in green to show the behavior of the 
visibilities and phases within individual CO lines. The positions 
of the CO lines are marked with the vertical ticks. The wave- 
length scale is in the heliocentric frame. A color version of this 
figure is available in the electronic edition. 




Fig. 4. Schematic view of our patchy model for Betelgeuse. a: 
Cross section of the model. The star is surrounded by two layers 
with an inhomogeneous velocity field, b: 3-D view of the model. 
The N- and E-axes define the plane of the sky. Only one layer is 
drawn for the sake of visual clarity. Within the cone specified by 
a half-opening angle of and the position (9, <p) with respect to 
the observer, the CO gas is assumed to be moving radially out- 
ward (or inward), while it is moving radially inward (or outward) 
in the remaining region. 
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Fig. 5. Comparison between our patchy model (T in = 2250 K, 2V 0Ut = 1 x 10 20 cm" 2 , v flow = 10 km s" 1 , = 60°, 6 = 40°, and 0=10°) 
and the AMBER data for Betelgeuse. In all panels, the solid lines represent the model, while the dots represent the observational 
data (data set #1). a: Normalized flux, b-d: Visibilities on the E0-G0-16m, G0-H0-32m, and E0-H0-48m baselines. The observed 
and model visibilities on the latter two baselines are binned with three and five pixels, respectively, e: Closures phase. The observed 
data and the model are binned with five pixels, f-h: Differential phases on the E0-G0-16m, G0-H0-32m, and E0-H0-48m baselines. 
The observed and model DPs on the latter two baselines are binned with three and five pixels, respectively. 
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Fig. 6. Images of Betelgeuse predicted by our patchy model 
shown in Fig. at different wavelengths in the CO line con- 
sisting of R(29) and R(12). a-d: Model images in the contin- 
uum, blue wing, line core, and red wing as indicated in the panel 
e. The extended MOLsphere is slightly visible in the panel c. 
e: Normalized flux (left ordinate) and visibility on the shortest 
baseline (right ordinate). The filled circles and solid line rep- 
resent the observed and model flux, respectively. The observed 
visibility and model prediction are shown by the filled diamonds 
and dashed line, respectively. The line positions are marked with 
the vertical ticks. A color version of this figure is available in the 
electronic edition. 
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Table 1. Summary of the AMBER observations of Betelgeuse and the calibrator Sirius. f b s : time of observation (Coordinated 
Universal Time=UTC), Tel.: telescope configuration, B p : projected baseline length, P. A.: position angle of the projected baseline on 
the sky, seeing in the visible, DIT: Detector Integration Time, Frames: the number of frames. 



# 


Name 


Night 


fobs Tel. B p P.A. 
(UTC) (m) O 


Seeing 

(") 


Airmass 


DIT 

(ms) 


Frames 


2008, HR.K (A/AA = 12000), 2.28-2.31 fim 


1 


Betelgeuse 


Jan. 08 


03:52:58 E0-G0-H0 16.00/31.97/47.96 73 


0.4 


1.19 


120 


2500 


Cal 


Sirius 


Jan. 08 


04:18:54 E0-G0-H0 16.00/31.99/47.99 72 


0.3 


1.01 


120 


2500 


2 


Betelgeuse 


Jan. 08 


04:48:24 E0-G0-H0 15.70/31.37/47.06 71 


0.3 


1.28 


120 


2500 


2006, MR_K (A/AA = 1500), 2.1-2.2 ^m 


3 


Betelgeuse 


Feb. 10 


02:45:08 E0-G0 15.45 70 


1.1 


1.33 


48 


6000 


Cal 


Sirius 


Feb. 10 


04:50:22 E0-G0 12.73 81 


1.1 


1.32 


48 


2500 
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Table 2. Summary of the VINCI observations of Betelgeuse. 
f b s : time of observation, Stations: siderostat stations, B v : pro- 
jected baseline length, P. A.: position angle of the projected base- 
line on the sky, and Visibilities: the result based on the wavelet 
analysis. 



Table 2. Continued 



# 


^obs 


Stations 


5 P 


P.A. 


Visibility 




(UTC) 




(m) 


o 




2001 Oct. 12 


1 


08:11:44 


E0-G0 


14.71 


74 


0.0988 + 0.0021 


2002 Dec. 06 


2 


05:10:46 


B3-C3 


7.70 


74 


0.4875 + 0.0066 


3 


05:37:38 


B3-C3 


7.89 


73 


0.4637 ± 0.0086 


2003 Jan. 07 


4 


02:40:04 


B3-C3 


7.50 


74 


0.4752 + 0.0082 


2003 Jan. 09 


5 


01:41:42 


B3-C3 


6.79 


74 


0.5653 ± 0.0069 


6 


02:10:29 


B3-C3 


7.25 


74 


0.5160 + 0.0063 


2003 Jan. 13 


7 


02:02:03 


B3-C3 


7.25 


74 


0.5128 + 0.0064 


8 


02:05:11 


B3-C3 


7.32 


74 


0.5093 + 0.0063 


9 


02:11:44 


B3-C3 


7.39 


74 


0.5060 + 0.0063 


10 


02:15:32 


B3-C3 


7.44 


74 


0.5060 + 0.0063 


11 


02:19:20 


B3-C3 


7.48 


74 


0.4989 + 0.0062 


12 


03:23:52 


B3-C3 


7.96 


73 


0.4503 + 0.0056 


13 


03:27:39 


B3-C3 


7.97 


73 


0.4461 +0.0055 


14 


03:31:28 


B3-C3 


7.97 


72 


0.4463 + 0.0055 


15 


04:00:31 


B3-C3 


7.97 


72 


0.4546 + 0.0057 


16 


04:04:18 


B3-C3 


7.96 


71 


0.4519 + 0.0057 


17 


04:08:07 


B3-C3 


7.95 


71 


0.4476 + 0.0056 


18 


04:33:46 


B3-C3 


7.83 


70 


0.4614 ± 0.0058 


19 


04:35:47 


B3-C3 


7.81 


70 


0.4666 ± 0.0058 


20 


04:39:36 


B3-C3 


7.78 


70 


0.4719 ±0.0059 


2003 Jan. 14 


21 


00:49:57 


B3-C3 


6.04 


73 


0.6780 ± 0.0068 


22 


00:54:02 


B3-C3 


6.13 


73 


0.6691 ±0.0068 


23 


00:58:04 


B3-C3 


6.21 


73 


0.6560 ± 0.0067 


24 


01:09:58 


B3-C3 


6.46 


73 


0.5993 ± 0.0055 


25 


01:14:00 


B3-C3 


6.53 


74 


0.5942 ± 0.0054 


26 


01:18:01 


B3-C3 


6.61 


74 


0.5870 ± 0.0054 


27 


01:27:23 


B3-C3 


6.78 


74 


0.5779 ± 0.0055 


28 


01:31:26 


B3-C3 


6.85 


74 


0.5802 ± 0.0057 


29 


01:41:48 


B3-C3 


7.02 


74 


0.5478 ±0.0051 


30 


01:45:52 


B3-C3 


7.08 


74 


0.5432 ±0.0051 


31 


01:49:54 


B3-C3 


7.14 


74 


0.5350 ± 0.0052 


32 


02:54:37 


B3-C3 


7.84 


73 


0.4871 ±0.0065 


33 


02:58:40 


B3-C3 


7.86 


73 


0.4811 ±0.0059 


34 


03:02:44 


B3-C3 


7.89 


73 


0.4940 ± 0.0057 


35 


03:14:34 


B3-C3 


7.94 


73 


0.4534 ± 0.0043 


36 


03:18:37 


B3-C3 


7.95 


73 


0.4534 ± 0.0043 


37 


03:22:39 


B3-C3 


7.96 


73 


0.4573 ± 0.0044 


38 


03:34:14 


B3-C3 


7.98 


72 


0.4649 ± 0.0045 


39 


03:38:18 


B3-C3 


7.99 


72 


0.4582 ± 0.0043 


40 


03:42:22 


B3-C3 


7.99 


72 


0.4550 ± 0.0044 


41 


04:20:35 


B3-C3 


7.88 


71 


0.4530 ± 0.0044 


42 


04:24:38 


B3-C3 


7.86 


70 


0.4643 ± 0.0043 


43 


04:28:43 


B3-C3 


7.83 


70 


0.4688 ± 0.0046 



# 


'obs 


Stations 


Bp 


P.A. 


Visibility 




(UTC) 




(m) 


o 




2003 Jan. 15 


44 


00:52:16 


B3-C3 


6.17 


73 


0.6375 ± 0.0087 


45 


00:56:20 


B3-C3 


6.26 


73 


0.6444 ± 0.0088 


46 


01:00:21 


B3-C3 


6.34 


73 


0.6269 ± 0.0086 


47 


01:14:12 


B3-C3 


6.61 


74 


0.6642 ± 0.0109 


48 


01:18:14 


B3-C3 


6.68 


74 


0.5902 ± 0.0080 


49 


01:31:23 


B3-C3 


6.92 


74 


0.5704 ± 0.0075 


50 


01:35:30 


B3-C3 


6.98 


74 


0.5612 ± 0.0074 


51 


01:45:14 


B3-C3 


7.15 


74 


0.5252 ± 0.0069 


52 


01:50:55 


B3-C3 


7.22 


74 


0.5228 ± 0.0068 


53 


01:56:38 


B3-C3 


7.30 


74 


0.51 10 ±0.0067 


54 


02:07:02 


B3-C3 


7.43 


74 


0.4975 ± 0.0065 


55 


02:12:43 


B3-C3 


7.51 


74 


0.4847 ± 0.0063 


56 


02:18:24 


B3-C3 


7.57 


74 


0.4898 ± 0.0064 


57 


03:57:19 


B3-C3 


7.96 


71 


0.4525 ± 0.0059 


58 


04:03:01 


B3-C3 


7.94 


71 


0.4555 ± 0.0059 


59 


04:08:42 


B3-C3 


7.91 


71 


0.4618 ±0.0060 


60 


04:19:11 


B3-C3 


7.87 


70 


0.4578 ± 0.0060 


61 


04:25:17 


B3-C3 


7.83 


70 


0.4765 ± 0.0063 


62 


04:30:58 


B3-C3 


7.78 


70 


0.4876 ± 0.0066 


2003 Jan. 18 


63 


01:11:11 


B3-C3 


6.76 


74 


0.5671 ±0.0106 


64 


01:14:20 


B3-C3 


6.86 


74 


0.5625 ± 0.0104 


65 


01:28:58 


B3-C3 


7.14 


74 


0.5296 ± 0.0097 


66 


01:48:44 


B3-C3 


7.41 


74 


0.5002 ± 0.0092 


67 


02:03:30 


B3-C3 


7.60 


74 


0.4676 ± 0.0086 


68 


03:13:00 


B3-C3 


7.97 


72 


0.4430 ±0.0081 


69 


03:15:42 


B3-C3 


7.98 


72 


0.4389 ±0.0081 


70 


03:52:06 


B3-C3 


7.94 


71 


0.4433 ± 0.0082 


71 


03:55:13 


B3-C3 


7.93 


71 


0.4458 ± 0.0083 


72 


04:00:32 


B3-C3 


7.89 


71 


0.4513 ±0.0082 


73 


04:12:01 


B3-C3 


7.81 


70 


0.4557 ± 0.0083 
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Table 3. Calibrators used for the VINCI observations of 
Betelgeuse. Name of the calibrators, spectral type (Sp.Type), 
uniform-disk diameters (duo), which are taken from Richichi & 
Percheron (2005 ). The number of the observations of the indi- 
vidual calibrators on each night is given in braces. 



Name 


Sp.Type 


rfuD (mas) 


Night 


e Peg 


K2Ib 


7.70 ± 0.24 


2001 Oct. 12 (4) 


/3Cet 


G9II-III 


5.18 ±0.05 


2001 Oct. 12 (6) 








2002 Dec. 06 (2) 


e Lep 


K4III 


5.90 ± 0.06 


2001 Oct. 12 (3) 








2002 Dec. 06 (1) 








2003 Jan. 09 (1) 


Sirius 


A1V 


5.60 ±0.15 


2001 Oct. 12 (4) 








2002 Dec. 06 (1) 








2003 Jan. 07 (1) 








2003 Jan. 18 (6) 


c5CMa 


F8Iab 


3.29 ± 0.46 


2001 Oct. 12 (2) 








2003 Jan. 07 (1) 


ix Eri 


M1III 


4.8 ±0.5 


2002 Dec. 06 (1) 


T Pup 


K1III 


4.38 ± 0.07 


2002 Dec. 06 (1) 








2003 Jan. 09 (1) 


Procyon 


F5IV-V 


5.37 ± 0.11 


2002 Dec. 06 (1) 








2003 Jan. 07 (1) 


1 Pup 


K5III 


3.8 ±0.4 


2002 Dec. 06 (1) 


^-Phe 


K5III 


2.69 ± 0.03 


2002 Dec. 06 (1) 


t Cet 


K1.5III 


3.28 ± 0.04 


2003 Jan. 07 (1) 


31 Ori 


K5III 


3.55 ± 0.06 


2003 Jan. 07 (1) 


ffCMa 


K4III 


4.13 ± 0.4 


2003 Jan. 07 (1) 


jSCnc 


K4III 


4.88 ± 0.03 


2003 Jan. 07 (1) 








2003 Jan. 09 (2) 


t Hya 


G9II-III 


3.1 ±0.2 


2003 Jan. 07 (1) 








2003 Jan. 09 (1) 


a Hya 


K3II-III 


9.44 ± 0.9 


2003 Jan. 07 (1) 








2003 Jan. 09 (6) 








2003 Jan. 13 (6) 


( Hya 




3.41 ± U.Uj 


2UU3 Jan. U/ (I) 








2003 Jan. 09 (4) 








2003 Jan. 13 (4) 


H Hya 


K2.5III 


4.69 ± 0.5 


2003 Jan. 09 (3) 








2003 Jan. 13 (9) 








2003 Jan. 15 (2) 


y Leo 


KHIIb 


7.70 ± 0.70 


2003 Jan. 09 (3) 


/?Ori 


B8Iab 


2.43 ± 0.05 


2003 Jan. 13 (16) 








2003 Jan. 14 (22) 








2003 Jan. 15 (15) 


V337 Car 


K3IIa 


5.09 ± 0.06 


2003 Jan. 14 (6) 


N Vel 


K5III 


6.92 ± 0.08 


2003 Jan. 07 (1) 



